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We propose an alternative mechanism for intracellular cargo transport which results from motor 
induced longitudinal fluctuations of cytoskeletal microtubules (MT). The longitudinal fluctuations 
combined with transient cargo binding to the MTs lead to long range transport even for cargos 
and vesicles having no molecular motors on them. The proposed transport mechanism, which we 
call "hitchhiking", provides a consistent explanation for the broadly observed yet still mysterious 
phenomenon of bidirectional transport along MTs. We show that cells exploiting the hitchhiking 
mechanism can effectively up- and down-regulate the transport of different vesicles by tuning their 
binding kinetics to characteristic MT oscillation frequencies. 
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Molecular motor mediated transport along micro- 
tubules (MTs) is a well studied phenomenon in vitro. 
There is an increasing number of studies of classic mi- 
crotubule motors like kinesin and dynein that shed light 
on their mechanochemistry in the idealized situation of 
in vitro single molecule assays [lj. Despite significant 
in vitro advances, understanding how intracellular trans- 
port works in vivo still remains one of the big challenges 
in molecular biology. Questions like how cellular cargo 
vesicles find their way through the cytoplasm and get 
targeted to their temporary or final destinations are at 
the heart of the problem. One of the major mysteries in 
this context is a phenomenon called "bidirectional trans- 
port" (BDT) The majority of cargos in the cell 
move in a bidirectional and remarkably symmetric man- 
ner. Despite the known kinetic and dynamic asymmetry 
of the underlying + and — end directed motors (kinesin 
and dynein respectively) , the vesicles seem to move with 
the same rates and run length distributions, and exhibit 
identical stalling forces, in each direction 0,0- 

Previous work on BDT has been dominated by the 
search for a hypothesized coupling element coordinat- 
ing the actions of motors attached to a cargo vesicle to 
produce the observed back and forth symmetric motion. 
But after years of quest for the molecules that establish 
the predicted coupling between dynein and kinesin, the 
molecular mechanism of the putative coordination model 
remains obscure. In its very basis the coordination 
model remains controversial, as it is hard to imagine any 
molecular regulator that would symmetrize the behavior 
of two so distinct motor species like dynein and kinesin. 
Nevertheless there is strong evidence from several model 
organisms that when either dynein or kinesin are 
impeded in their function the bidirectional transport gets 
affected in both + and — directions symmetrically. Here 
we propose a physically and biologically plausible molec- 
ular model for bidirectional transport that can account 
for the symmetric behavior in both directions. 

Our model has two basic ingredients: 1. Motors of 
either or both types (not necessarily localized to 



the cargo) generate longitudinal MT fluctuations 

by already known mechanisms, e.g. the one depicted in 
Fig. la. This motor induced confined MT sliding we call 
"jabberwalking" of the underlying motors. While inter- 
MT shearing forces are best known for motile organelles 
containing the axoneme structure Q , an increasing num- 
ber of studies show vigorous sliding of cytoplasmic MTs 
when they become detached from cellular structures 0. 
Two-point microrheology studies of cytoskeletal stress- 
strain fluctuations that demonstrate the presence active 
force doublets go along the same line of evidence. 
2. Vesicles transiently couple to the random, 
stochastically oscillating MTs and attain their speed 
temporarily. As a consequence, the vesicles experience 
strongly enhanced diffusion. This process we call "hitch- 
hiking" of the vesicles. The occurrence of such a mecha- 
nism is most clearly observed for reticulopodial cytoskele- 
ton extensions and chlamydomonas flagella 0. In 
those systems, artificial and endogenous cargos (Fig. lc) 
are found to move as a consequence of motor driven MT 
sliding in a bidirectional manner. More recently, evidence 
was found for hitchhiking in S2 drosophila cell lines 
indicating a possible ubiquity of the process. 

Below we first focus on the hitchhiking process and 
compute how a vesicle transiently coupling to a single 
longitudinally oscillating MT moves on long timescales 
(Fig. lb). Each time the cargo binds to the MT, it as- 
sumes the latter's direction and velocity v — wmt! when it 
is unbound, we assume complete immobility i.e., v = 0. 
We neglect here the thermal diffusion of the unbound 
vesicle due to its smallness compared to the active trans- 
port we consider here (see discussion below). The veloc- 
ity vmt (t) itself is a random process with the particular 
property that xmt (t) = J Q vmt (t) dr is a bounded vari- 
able, i.e., we assume the MT to move in spatially con- 
fined fashion. Consequently a cargo simply permanently 
sticking to the MT does not get far, and so moves in a 
confined manner as well. However by virtue of the switch- 
ing (attachment/detachment) process, the cargo coordi- 
nate x (t) can become unbounded! 
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FIG. 1: a. Longitudinal MT noise caused by "jabberwalking" 
of double attached motors on elastically tethered MTs. 6. The 
hitchhiking mechanism: Vesicle binding and unbinding kinet- 
ics, combined with local MT oscillations, induces diffusive, 
long range transport. MT surface polarity can also give rise 
to a directed drift, c. Reticulomyxa transports microspheres 
(M) and organelles (O) via MT sliding (adapted from Q). 

The vesicle binding process is described by a 2-state 
random variable B (t) € {0, 1} with the characteristics of 
telegraphic noise : B (t) equals 1 if the cargo sticks 
to the MT or for a detached resting cargo. We de- 
note the stochastic switching rates between those states 
by k oS and fc on . We next write the velocity of the 
cargo vesicle as a composite random process v (t) = 
B (t) vmt (t) ■ The cargo velocity correlation function is 
given by (v (ii) v (t 2 )) = (B (h) B (t 2 ) v M t "mt (t 2 )), 
where (• • •) denotes the ensemble average over all real- 
izations of B and v. For simplicity, at first we focus on 
the limiting case of velocity independent rates k and 
k on (we consider the general case of velocity dependence 
below). In this case, the vesicle binding becomes sta- 
tistically independent from the MT motion (Bvmt) = 
(B)(vmt), B (t) becomes a standard (asymmetric) 
Markovian telegraphic noise process [ll], and we easily 
compute (v(ti)v(t 2 )) = Cmt (h, t 2 ) Cb (h — t 2 ) with 
Cmt(*i,*2) = (umt(*i)«mt(*2)) and C B {h-t 2 ) = 
K 2 +K(1-K) e -2|ti-* 2 |/* B w_ Tne tw0 constants K and 
t sw characterize the binding behavior. K = (B) = 
k on / (fc on -(- fc off ) denotes the equilibrium binding con- 
stant and t sw — 2/ (fc on + k ) is the "mean switching 
time." Note that whereas the switching process B (t) is 
assumed to be stationary in the statistical sense, the ve- 
locity vmt (t) can in general be a non-stationary stochas- 
tic or even a purely deterministic process, and its auto- 
correlation function Cmt (ti,t 2 ) is not necessarily time 



homogeneous. The mean-square-displacement (msd) of 
such a vesicle is obtained from 

((x (t) - x (Q)) 2 ) =J J CMT(tuh)C B (t 1 -t 2 )dt 1 dt2 

(1) 

This general expression relates the msd of the cargo to 
the fluctuations of a single microtubule and the bind- 
ing/unbinding kinetics of the cargo to the MT. For the 
special case of a statistically stationary process v% (t) , the 
correlator Cmt (ti,t 2 ) = Cmt [ti — t 2 ) becomes homoge- 
neous in time and we can simplify Eq.^ further 

((x(i)-z(O)) 2 ) = 2^ (t - t) C b (r) C MT (r) dr (2) 

Now we illustrate our general formulas by focusing on two 
particular possible MT shaking processes: (1) Determin- 
istic oscillation of the MT with vmt (i) — V sin cut, and 
(2) An idealization of MT motion as an overdamped har- 
monic oscillator driven by Markovian telegraphic noise 
(representing the action of motors, Fig. la). 

Periodically oscillating MT. In the first case, we have 
Cmt (ii; t 2 ) = V 2 sin iot\ sin ut 2 , and it is easy to evalu- 
ate the rhs of Eq.Q] (Eq.|5| cannot be used because of the 
non-stationarity of v\ (t)). The result can be simplified 
in two limiting cases. On short timescales t <ti t sw , the 
vesicle does not have enough time to bind/unbind from 
the MT, and we simply have Cb (t) ~ K. From Eq.Q] 

we obtain /(x (t) — x (0)) 2 ^} = Klu~ 2 V 2 (1 — costw) 2 , as 

expected for a particle strictly following the MT with 
probability K. In particular, the vesicle stays spatially 
confined. In the opposite limit t 3> t sw the vesicle mo- 
tion becomes diffusive: l(x (t) — x (O)) 2 ^) = const. + 2Dt, 
with the diffusion constant given by 

D = K (1 - K) V 2 t sw (4 + coh^y 1 (3) 

This remarkable expression says that a vesicle stochasti- 
cally coupling to a periodically oscillating MT diffuses 
with an efficiency that depends on the fine tuning of 
the MT oscillation frequency w and the stochastic vesicle 
switching time £ sw . The long time (t 3> t sw ) transport 
efficiency is maximized for K opt — 1/2 and t° pt = 2lu~ 1 1 
and falls off to zero away from these values. The intu- 
itive meaning of the first result (K° pt = 1/2) is clear: 
If the vesicle sticks too strongly to the MT (K = 1) 
or not at all (K = 0), there is no long range transport 
as it either moves with the MT (in a confined manner) 
or not at all. The optimum occurs for an intermediate 

value. The second finding t° p ^ = 2uj~ 1 , which resembles 

° sw I — I ' 

stochastic resonance phenomena [lfj, also has a simple 
interpretation: If the vesicle takes too long a ride on the 
MT (t S w w_1 ), hs average displacement cancels be- 
cause of the pure "back and forth" motion of the MT. 
If the ride is too short (t sw <C uj _1 ), a similar argument 
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applies. In this rather trivial example, we already see an 
interesting theme: Optimal transport of vesicles requires 
a fine tuning of vesicle binding and MT oscillations. 

MT driven by stochastic motor noise. In a second more 
realistic approach, we model longitudinal MT oscillations 
by an overdamped Langevin equation with a harmonic 
restoring force coming from MT attachment /confinement 
(Fig. la,b): £i M T = ~Cx M t + F mot (±uT,t) ■ Here £ is 
the MT longitudinal friction constant and C the MT 
restoring spring constant. The actively generated mo- 
tor force -F m ot (%t , t) depends on the detailed motor 
mechanochemistry, which we effectively model by a lin- 
ear force-velocity relation v mot (Fmot) /vq = 1 — F mo t/F 
with two parameters: Vo, the maximal (zero-load) ve- 
locity, and Fq, the motor stalling force. Combining this 
with the Langevin equation and ±mt = Vmot yields the 
effective equation of MT motion: 

Coff (f) ±MT = -CXMT + Fq (t) (4) 

£eff (t) = £, + Fq (t) /vq (t) is the effective friction constant. 

EqJH states that the motors contribute to an effec- 
tive external force Fq (t), but also give rise to increased 
effective friction £ e g > £. The driving motor force 
Fq (t) and velocity vq (t) are stochastic variables, which 
can switch between two values. The dynamics of this 
switching generally depends on the MT-motor attach- 
ment geometry. In the simplest arrangement (Fig. la), 
motors bind rigidly and run actively on both MTs in a 
symmetric manner. In this case Fq (t) and vq (t) both 
switch between two values, which for simplicity we as- 
sume to be equal in magnitude but of opposite sign, i.e.. 
Fq (t) — ±Fo and vq (t) — ±vq (same number of motors 
of same strength on both sides), which results in a time 
independent friction constant £ e ff — £ + Fq/vq. We as- 
sume that the motors stochastically switch direction with 
an exponentially distributed switching time, i.e., Fq (t) 
is described by symmetric Markovian telegraph noise 
[HI with p {F Q (t) = ±F ) = 1/2 and (F (h) F Q (i 2 )) = 
Fq exp (—2 \ti — t2 1 /T p ). Here T p is the processivity 
time of the motors (average time between direction 
changes) . 

In long time limit, where imt in Eq.0] becomes a 
stationary process, we can exploit Eq.|21 to evaluate 
the msd of the vesicle, provided that we can compute 
Cmt (t) = (xmt (t + t) xmt (£))• To accomplish this, 
we use the solution of Eq.|U in the limit i~^Ct ^> 1: 
x(t) = Jo Cff F o(rle-^ lc(t - T) dr. After some calcula- 
tion, this leads to |l2j: 

The two characteristic timescales are now the MT relax- 
ation time T r = £ c ff/C and the motor processivity time 
T p . Inserting Eq.0 into Eq.H gives a lengthy expression 



which in the limit t max(T p ,T r ,t sw ) can be used to 
compute the long-time diffusion constant 

D = Fq 2 2K (1 - K) t sw T p 

{2T r ~t- i S w) (^sw "I - Fp) (Tp + 2F r ) 

As in the previous example (periodic MT oscillation), 
the vesicle mobility is always maximized for K° pt = 1/2, 
i.e., for an intermediate binding strength. The opti- 
mal vesicle switching time t°P* = ^J2T T T V is also eas- 
ily obtained from Eq.[6| Interestingly, for fixed T r , Eq.[B] 
predicts the optimal ratio T p /T r = 2. For those val- 
ues we obtain the maximal diffusion constant -D max = 
Fq (£ + Fq/vq) 1 C~ 1 /32, i.e., stronger and faster mo- 
tors {Fq and vq large), weaker MT confinement (C small) 
and smaller MT friction give rise to more efficient trans- 
port. 

Multiple MTs. For n different oscillating MTs, the 
particle velocity is v(t) = YZ=i 5 k,B(t)VMT k (t) with 
VMT k the velocity of the fcth MT, 6k.i the Kroneker-delta 
and B (t) € {0, 1, 2, .., n} the "binding variable" which in- 
dicates to which MT the vesicle is bound at time t (0 rep- 
resents the unbound state). For the two-MT situation in 
Fig. la, n = 2 and VMT t = —vmt 2 (anticorrelated MT ve- 
locities). Using Eq.JSJas before yields the resonance con- 
dition: fc on > fc off , k oS = V2{T r T p y 1/2 and T p = 2T r , 
which differs slightly from the single MT case because it 
is now more favorable to jump between the two MTs than 
to spend time in the unbound state. The corresponding 
diffusion constant D max = Fq (£ + Fq/vq)^ 1 C~ 1 /8 is 4 
times larger than in the single MT case. This indicates 
that the efficiency of transport, as well as the natural 
strategy for optimizing it, can depend on the effective 
number of participating MTs. 

Hitchhiking vs thermal diffusion. In the limit of motor 
forces larger than viscous forces (F » £,v ), we obtain 
the rough estimate D max w \{Fq/C)vq. Typical MT 
sliding velocities vq w 1-5 /im/s and oscillation ampli- 
tude Fq/C w 0.1-1/im yield F» max w 0.01-0.5Mm 2 s _1 . In 
comparison, a typical organelle with diameter 500nm ex- 
periences a large cellular viscosity r\ ~ 0.3 Pa • s |l( and 
has a thermal diffusion constant Dth » 3 • 10~ 3 /im 2 s _1 . 
On the other hand, nanometer sized molecules experi- 
ence a much smaller effective viscosity (close to that of 
water r\ w 10~ 3 Pa-s) and so have D t h « 10-100 /j,m 2 s _1 . 
Therefore large objects like vesicles and organelles can 
strongly benefit from hitchhiking while smaller molecules 
are more efficiently transported by thermal diffusion. 

Hydrodynamic stress and biased hitchhiking. It is 
straightforward to show that the hitchhiking process as 
described above for large t becomes an unbiased diffusive 
process even for an asymmetric MT shaking . We con- 
sider next an interesting generalization of the hitchhiking 
model which gives rise to biased directed transport. 

In the following, we drop the previous assumption 
(Bvmt) = (B) (vmt) i-e., the statistical independence 



4 



of the MT velocity and the vesicle-MT binding pro- 
cess. For example, a statistical coupling of binding and 
MT sliding can appear when the off rate k oS (vmt) be- 
comes velocity dependent due to a hydrodynamic drag 
force (with respect to the essentially stationary cyto- 
plasm) acting on a moving vesicle. The presence of this 
force can break the symmetry of transport in two dif- 
ferent ways. The general form of the off rate will be 
k off (v M t) /k oS (0) w l + C2 v MT + c 3 v MT + - ■ ■ . The linear 
term must vanish, if we make the plausible requirement 
that k («mt) is minimal for vmt = (no hydrodynamic 
stress = maximal binding strength) . We assume that the 
two remaining terms c 2 v MT and C3« MT are small. We call 
the coefficient c 2 the dynamical bias coefficient, be- 
cause it gives rise to a particle drift only if the MT has 
different forward and backward velocity. We call C3 the 
polarity bias, as it can give rise to drift even for a time 
reversal symmetric dmt (*)• Physically, this term stems 
from the polarity of MT surface and the resulting polar- 
ity of the interaction with the vesicle. 

We next derive the mean drift velocity (v (t)) of 
the vesicle in the simplifying limiting case of rapid 
binding equilibration i.e., k on + k oS ^> ^MT.char 
where WMT,char denotes the characteristic oscillation fre- 
quency of the MT. In this limit, the vesicle bind- 
ing state equilibrates at each instant of time and 
the conditional binding probability becomes a function 
of the instantaneous MT velocity: p(B — 1\vmt) — 
*°7 (fc off (vmt) + fc on ) « Ko ~ K (1 - K ) (02i& T + 
c 3 v MT ) with K = fc°7 (fc off (0) + fc on ) . Conse- 
quently the vesicle's mean drift velocity simplifies to 
(v(t)} = («mt(*)p(£ = 1|umt(*))) = (vmt (t))K - 
K (1 - Ko) (c 2 (v MT (t))+c 3 (u MT (t)) and its long time 
mean (vesicle drift) reads: 

v = -Ko (1 - Ko) (c 2 (<tT + c 3 « T )) (7) 

where the long time average / = limr— >oo T~ l J Q T / (t) dt. 
{vmt) vanishes because of MT confinement. 

As an example, we consider a (deterministic) asym- 
metric square wave MT oscillation with vmt (t) period- 
ically switching between vmt — V\ for a time T\ and 
vmt = V 2 for time T 2 with V\ V 2 < and T x Vi +T 2 V% = 
(zero mean). For such a choice from Eq.|7| we ob- 
tain v = -K (1 - K ) (Tj + T 2 y x (c 2 (T X V? + T 2 V 2 3 ) + 
C3 (T\Vi + T 2 V 2 )) which in general has both the dynam- 
ical (c 2 ...) and the polar (03...) drift contributions. For a 
completely symmetric shaking (V 2 = —V\, T\ = T 2 ), the 
dynamical part vanishes as expected, whereas the polar 
part always stays present with a sign opposite to c 3 . 

In conclusion, we have outlined a transport mechanism 
alternative to the standard "cargo hauled by a motor" 
model. Active longitudinal motions of MTs, combined 
with even weak non-specific cargo— MT binding, natu- 
rally lead to this transport scenario. Remarkably, even 
small scale (spatially confined) MT oscillations induce 



a long distance transport on the cellular scale. Hitch- 
hiking on approximately symmetric MT bundle arrange- 
ments as in Fig. la inherently bears the feature of veloc- 
ity and run-length distribution symmetry often observed 
in BDT 0, Q. Very large moment velocities US 
(up to > 10 /im/s) — another characteristic of BDT — may 
be attributable to the enhanced motor efficiency known 
for cooperatively synchronized motor batteries perform- 
ing filament sliding 13] and elastic MT relaxations 
Other signature behaviors of our mechanism, such as vis- 
ible sliding of MTs, tandem motion of vesicles, exponen- 
tial velocity relaxations and strong MT bending deforma- 
tions have been observed recently 0, El • In general, we 
expect a combination of hitchhiking and hauling mecha- 
nisms to be present in varying proportions for different 
cargos and cells. Future experiments might reveal the rel- 
ative importance of their contributions and the "master 
plan" behind them. 
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